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Abstract 

The multisite phosphorylation-dephosphorylation cycle is a motif repeatedly used in cell signaling. 
This motif itself can generate a variety of dynamic behaviors like bistability and ultrasensitivity without 
direct positive feedbacks. In this paper, we study the number of positive steady states of a general 
multisite phosphorylation-dcphosphorylation cycle, and how the number of positive steady states varies 
by changing the biological parameters. We show analytically that (1) for some parameter ranges, there 
are at least n + 1 (if n is even) or n (if n is odd) steady states; (2) there never are more than 2n — 1 
steady states (in particular, this implies that for n = 2, including single levels of MAPK cascades, there 
are at most three steady states); (3) for parameters near the standard Michaelis-Menten quasi-steady 
state conditions, there are at most n + 1 steady states; and (4) for parameters far from the standard 
Michaelis-Menten quasi-steady state conditions, there is at most one steady state. 

Keywords: futile cycles, bistability, signaling pathways, biomolecular networks, steady states 

1 Introduction 

A promising approach to handling the complexity of cell signaling pathways is to decompose pathways into 
small motifs, and analyze the individual motifs. One particular motif that has attracted much attention in 
recent years is the cycle formed by two or more inter-convertible forms of one protein. The protein, denoted 
here by So, is ultimately converted into a product, denoted here by S n , through a cascade of "activation" 
reactions triggered or facilitated by an enzyme E; conversely, S n is transformed back (or "deactivated") 
into the original So, helped on by the action of a second enzyme F. See Figured! 





Figure 1: A futile cycle of size n. 



Such structures, often called "futile cycles" (also called substrate cycles, enzymatic cycles, or enzymatic 
inter-conversions, see [1]), serve as basic blocks in cellular signaling pathways and have pivotal impact on 
the signaling dynamics. Futile cycles underlie signaling processes such as GTPase cycles [2], bacterial 
two-component systems and phosphorelays [3,4] actin treadmilling [5]), and glucose mobilization [6], as 
well as metabolic control [7] and cell division and apoptosis [8] and cell-cycle checkpoint control [9]. One 
very important instance is that of Mitogen- Activated Protein Kinase ( "MAPK" ) cascades, which regulate 
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primary cellular activities such as proliferation, differentiation, and apoptosis [10-13] in eukaryotes from 
yeast to humans. 

MAPK cascades usually consist of three tiers of similar structures with multiple feedbacks [14-16]. 
Each individual level of the MAPK cascades is a futile cycle as depicted in Figure Q] with n = 2. Markevich 
et aVs paper [17] was the first to demonstrate the possibility of multistationarity at a single cascade 
level, and motivated the need for analytical studies of the number of steady states. Conradi et al. studied 
the existence of multistationarity in their paper [19], employing algorithms based on Feinberg's chemical 
reaction network theory (CRNT). (For more details on CRNT, see [31,32].) The CRNT algorithm confirms 
multistationarity in a single level of MAPK cascades, and provides a set of kinetic constants which can 
give rise to multistationarity. However, the CRNT algorithm only tests for the existence of multiple steady 
states, and does not provide information regarding the precise number of steady states. 

In [18], Gunawardena proposed a novel approach to the study of steady states of futile cycles. His 
approach, which was focused in the question of determining the proportion of maximally phosphorylated 
substrate, was developed under the simplifying quasi-steady state assumption that substrate is in excess. 
Nonetheless, our study of multistationarity uses in a key manner the basic formalism in [18], even for the 
case when substrate is not in excess. 

In Section O we state our basic assumptions regarding the model. The basic formalism and background 
for the approach is provided in Section [3j The main focus of this paper is on Section 01 where we derive 
various bounds on the number of steady states of futile cycles of size n. The first result is a the lower 
bound for the number of steady states. Currently available results on lower bounds, as in [29], can only 
handle the case when quasi-steady state assumptions are valid; we substantially extend these results to the 
fully general case by means of a perturbation argument which allows one to get around these restricted 
assumptions. Another novel feature of our results in this paper is the derivation of an upper bound 
of 2n — 1, valid for all kinetic constants. Models in molecular cell biology are characterized by a high 
degree of uncertainly in parameters, hence such results valid over the entire parameter space are of special 
significance. However, when more information of the parameters are available, sharper upper bounds can 
obtained, see Theorems H] and [5j We finally conclude our paper in Section [5] with a conjecture of an n + 1 
upper bound. 

We remark that the results given here complement our work dealing with the dynamical behavior 
of futile cycles. For the case n = 2, [25] showed that the model exhibits generic convergence to steady 
states but no more complicated behavior, at least within restricted parameter ranges, while [27] showed a 
persistence property (no species tends to be eliminated) for any possible parameter values. These papers 
did not address the question of estimating the number of steady states. (An exception is the case n = 1, 
for which uniqueness of steady states can be proved in several ways, and for which global convergence to 
these unique equilibria holds [27].) 

2 Model assumptions 

Before presenting mathematical details, let us first discuss the basic biochemical assumptions that go into 
the model. In general, phosphorylation and dephosphorylation can follow either distributive or processive 
mechanism. In the processive mechanism, the kinase (phosphatase) facilitates two or more phosphorylations 
(dephosphorylations) before the final product is released, whereas in the distributive mechanism, the kinase 
(phosphatase) facilitates at most one phosphorylation (dephosphorylation) in each molecular encounter. 
In the case of n = 2, a futile cycle that follows the processive mechanism can be represented by reactions 
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as follows: 



S + E < — ► ES < — ► ESi — ► S 2 + E 
S 2 + F< — > FS 2 < — > FSi — ► So + F; 

and the distributive mechanism can be represented by reactions: 

S + E < — * ESo — > Si + E < — ► ESi — ► S 2 + E 
S 2 + F < — ► FS 2 — > Si + F < — > FSi — ► So + F. 

Biological experiments have demonstrated that both dual phosphorylation and dephosphorylation in MAPK 
are distributive, see [14-16]. In their paper [19], Conradi et al. showed mathematically that if either phos- 
phorylation or dephosphorylation follows a processive mechanism, the steady state will be unique, which, 
it is argued in [19], contradicts experimental observations. So, to get more interesting results, we assume 
that both phosphorylations and dephosphorylations in the futile cycles follow the distributive mechanism. 

Our structure of futile cycles in Figure Q] also implicitly assumes a sequential instead of a random 
mechanism. By a sequential mechanism, we mean that the kinase phosphorylates the substrates in a 
specific order, and the phosphatase works in the reversed order. This assumption dramatically reduces the 
number of different phospho-forms and simplifies our analysis. In a special case when the kinetic constants 
of each phosphorylation are the same and the kinetic constants of each dephosphorylation are the same, 
the random mechanism can be easily included in the sequential case. Biologically, there are systems, for 
instance the auto-phosphorylation of FGF-receptor-1, that have been experimentally shown to follow a 
sequential mechanism [33]. 

To model the reactions, we assume mass action kinetics, which is standard in mathematical modeling 
of molecular events in biology. 



3 Mathematical formalism 

In this section, we set up a mathematical framework for studying the steady states of futile cycles. Let us 
first write down all the elementary chemical reactions in Figure [TJ 

k on fc t 
So + E^ESo ^°S! + E 

k oS 



&on n _i fc , 
S n -i + E < — ESo S n + E 
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^on i , 
S 1 + F7^FS 1 S + F 
l oS Q 



S n + F < — FS n — > S n -i + F 

l 0Sn-l 



where k O u , etc., are kinetic parameters for binding and unbinding, ESo denotes the complex consisting of 
the enzyme E and the substrate So, and so forth. These reactions can be modeled by 3n + 3 differential- 
algebraic equations according to mass action kinetics: 



ds 

dt 

dsj 

dt 

dCj_ 

dt 
ddk_ 
dt 



-k OIlo s e + /e offo c + i C at d i 

-kon.Sie + k Q gCi + k^. ^-i - lon^Sif + ^ff.^dj + katidi+i, i = 1, ■ ■ ■ ,n - 1 
konj Sje - (k oS . + k cat . ) Cj , j = 0, . . . , n - 1 (1) 

fon fc _i«k/ - (^offfc_i + Z cat fc _i)4, fc = 1, 
together with the algebraic "conservation equations" : 



n-l 



E tot = e + J2 Ci > 


n 

*tot = / + E^ ( 2 ) 
l 

n n—l n 

s tot = J2 Si+ J2° i+ J2 di - 



1 

The variables so, . . . , s n , Co, . . . , c n _i, di, . . . , d n , e, / stand for the concentrations of 

Sq, ... , S n , ESo, ■■■ , ES n -i,FS±, . . . , FS n , E, F 
respectively. For each positive vector 

k =(^0110, • • • > ^on„_i, ^off ' • • • ' ^off„_i' ^cat ' • • • ' ^cat n _!i 

lon ,- • • , U„_i, l o ff , ■ ■ ■ , J ff n -i' l cat > ■ ■ ■ > Z cat n _i) G K^ -6 



(of "kinetic constants") and each positive triple C = (-Etot' ^tot' ^tot)' we nave a different system C). 
Let us write the coordinates of a vector x G IR^ n+3 as: 

x = (sq, • • • , >s n , cq, • • • , c n _i, di, . . . , d n , e, /), 



and define a mapping 



$ : M 3n+3 x IR^" 6 xl 3 + — ► M 3n+3 
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with components $1, . . . , &3 n +3 where the first 3n components are 

*i(aj, k,C) = -A; ono soe + ^ ff Co + Z cato d i 
and so forth, listing the right hand sides of the equations ([1]), $3 n +i is 

n-l 





and similarly for and we use the remaining equations in (J2J). 

For each k, C, let us define a set 

2(/c,C) = {x|*(a;,As,C) = 0}. 

Observe that, by definition, given x € M^™ +3 , x is a positive steady state of S(re, C) if and only if x € C). 
So, the mathematical statement of the central problem in this paper is to count the number of elements 
in Z(k,C). Our analysis will be greatly simplified by a preprocessing. Let us introduce a function 

* : M 3n+3 x R+ n ~ 6 xR} — ► M 3n+3 

with components . . . , \&3 n +3 defined as 

= + $n+l 
*i = $i + $„ +J + $2n+i-l + * = 2, . . . ,n 

*j ; = , j = n + 1, . . . , 3n + 3. 

It is easy to see that 

Z{k,C) = {x\V(x,K,C) = 0}, 

but now the first 3n equations are: 

*i = k&U-A ~ ^cati-iQ-i = °> * = 1, ,n, 
tfWi+j = feon^Sj-e - (fc off . + fc cat ,)c,- = 0, j = 0, . . . , n - 1 

*2n+fc = lon^Skf ~ ('off,-! + l c&i k ^) d k = 0, fc = 1, . . . ,n, 
and can be easily solved as: 

Sj+i = Xi{e/f)si, (3) 

(4) 

''+■!. ~ — • ( 5 ) 



K Mi 

fSi+l 



where 



. k ca,ti L Mi k caU + A; offi Z cat + Z off 
Ai = -rr-^ , K Mi = ; > ^M, = ; ; » = 0, . . . , n - 1. (6) 
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We may now express Sj, 1 c * an< ^ Si di in terms of sq, k, e and /: 



^ s i = s o f 1 + A o (y) + A oAi (y) + • • • + A • • • A n _i ^ j 

?— •fe + ^(j) + - +i ^ ! (7)")— »(?)■ 



?*-'-fe(7) + ^(7)" + - + ^(7)>'-*(7 

Although the equation \& = represents 3n + 3 equations with 3n + 3 unknowns, next we will show that 
it can be reduced to two equations with two unknowns, which have the same number of positive solutions 
as \P = 0. Let us first define a set 

S(k,C) = {(u,v) G M+ x R + \Gi' e (u,v) = 0,G%' C (u,v) = 0}, 

where G* ,C , G*^ C : R 2 ^ — > R are given by 



G?' C («,«) = « («^f («) - ^(^BtotMot) " ^tot/^tot + 
G£' C (u,i;) = <f$(u)<p$(uy + (<^(u) - S tot ^(u) + F tot u^(u) + F t ot v - 5 tot . 

The precise statement is as follows: 

Lemma 1 There exists a mapping 5 : R 3n+3 — ► R 2 such that, for each k,C, the map (5 restricted to 
Z(n,C) is a bijection between the sets Z(n,C) and S(k, C). 

Proof. Let us define the mapping 5 : R 3n+3 — ► R 2 as 5(x) = (e/f, sq), where 

x = (s ,...,s n ,c ,...,c n _i,di,...,ci n ,e, /). 

If we can show that 5 induces a bijection between Z(k,C) and S(k,C), we are done. 

First, we claim that 5(Z(k,C)) C <S(k, C). Pick any x G Z(k,C), we have that x satisfies flS}-©. 
Moreover, $3 ra+ 2(x, «, C) =0 yields 

E tot = e + es iPi(j), 

and thus 

„ _ ^tot 



l + s oV f{e/f) 

Using $ 3n+ i(x,K,C) = and $ 3n+2 (x,K,C) = 0, we get: 

gtgt _ e(l + S0 ^(e//)) 
F tot f(l + 8 <p5(e/f)y 

which is G*' (e/f, sq) = after multiplying by 1 + so(p%(e/ /) and rearranging terms. 
To check that G%'(e/f, sq) = 0, we start with $3 n+ 3(x, k,C) = 0, i.e. 



(8) 



0) 



ra-1 

s tot = Yl Si + Yl Ci + X! ^- 



1 
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Using ([7|) and (|HJ), this expression becomes 

Q - ~ ^r 6 ^ I gtotfOyjWj I F tot S 0^2( e //) 

Stot " ^o( 7 ) + 1 + S0 ^ (e//) + 1 + s^e/f) 

* ( e \ egtotgoPjKf//) *tot s o¥>2( e //) 
/(l + «o^(e//)) 1 + so^(e//) ' 

where the last equality comes from Q. 

After multiplying by 1 + SQip%(e/ /), and simplifying, we get 



¥>o(y)v>2 (y)*o + ( Po( j) " Stottfij) + j^otVi(j) + ^tot^Wj *o - 5 tot = 0, 

that is, G2 ,C (e/f, sq) = 0. since both G\' C (e/ f , sq) and G% (e/ f, sq) are zero, <5(x) G cS(k,C). 

Next, we will show that S(k,C) C <5(£(k, C)). For any y = (u, v) £ <S(k,C), let the coordinates of x be 



defined as: 



s = v 
Si+i = XiUSi 

e = Etot 



1 + s cp1(u) 

f= e - 

u 

esi 

Ci = 



<k+l — 



K Mi 

fSi+l 



L 



Mi 



for i = 0, . . . , n — 1. It is easy to see that the vector x = (so, . . . , s n , Co, . . . , c n _i, d\, . . . , d n , e, /) satisfies 
<3?i(x, ft, C) = 0, . . . , $3 n+ i(x, k, C) = 0. If <I>3 n _|_2(x, k, C) and &3 n+ 3(x, k,C) are also zero, then x is an 
element of Z(k,C) with <5(a;) = y. Given the condition that G^' C (u, v) = (i = 1,2) and u = e/f,v = so, 
we have G^' (e//, sq) = 0) an d therefore Q holds. Since 



^tot 



" l + so<pt(e/f) 

in our construction, we have 

n 

F tot = f(l + So ^(e/f)) = f + J2di 

i 

To check <E>3 n+ 3(x, re,C) = 0, we use 

Gf(e/f,s ) = 



1 + s ^(e/f) 

as G£ ,c (e//,s ) = and 1 + s Q ip%(e/f) > 0. Applying we have 



V> e , ^V"^ - . ■ ~ K ( r i f~\ i eF tot^?(e//) jtotfg^(g//) _ o 
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It remains for us to show that the map 8 is one to one on Z(k,C). Suppose that 5(x 1 ) = 5(x 2 ) = (u,v), 
where 

x = (sq, . . . , s n , Cq, . . . , c ra _ 1 , d^, . . . , d n , e , / ), i = 1,2. 

By the definition of 5, we know that Sq = Sg and e 1 // 1 = e 2 // 2 . Therefore, sj = sf for i = 0, ... ,n. 
Equation flSJ) gives 

e i = gtpt = e 2 
l + v(fl(u) 

Thus, / 1 = / 2 , and cj = c 2 , = d 2 +1 for i = 0, . . . , n — 1 because of ©-(IS]). Therefore, x 1 = x 2 , and 8 
is one to one. I 

The above lemma ensures that the two sets Z(k, C) and S(k,C) have the same number of elements. From 
now on, we will focus on S(k,C), the set of positive solutions of equations G*' (u, v) = 0,G2 ,C (u,v) = 0, 
i.e. 

Gf(u, v) = v (u<p$(u) - rt(u)E tot /F tot ) - E tot /F tQt +u = 7 (10) 
Gf(u, v) = ^(u)^(u)v 2 + (<tf(u) - S toW %(u) + F tot utf(u) + F tot ^(u)) v - 5 tot = 0. (11) 



4 Number of positive steady states 

4.1 Lower bound on the number of positive steady states 



One approach to solving (fT0j) - (fTT|) is to view G^ (u, v) as a quadratic polynomial in v. Since 
equation (fTT]) has a unique positive root, namely 

„. _ -H^ c {u) + x /ff^(n) 2 + 45 tot V? g(n)^(n) ^ 



2^{u)^{u) 
where 

H K > c (u) = y$(u) - S toW %(u) + F to t«^(«) + ^totV2(«). (13) 
Substituting this expression for v into (fTUj) , and multiplying by </?q (it) , we get 



TTK,Cf..\ ._ V ; V V ; tOt^tn /^2V ^ / \ -^tot ,„K(„,\ \ f^t,.Kf.,\_ L .., n Kf,.\ _ n 

b \ u ) ■- Kl \ \ u ^i\ u ) - — ftW hp — VoW+^VoW - u - 

2( P2K U ) \ b tot J -^tot 

(14) 

So, any («, v) G S(k,C) should satisfy (fl~2l) and {H]). On the other hand, any positive solution u of (fT4"P 
(notice that ^Pq(u) > 0) and v given by (fT2|) (always positive) provide a positive a solution of (rf0j) - (fTT]) . 
that is, (u,v) is an element in S(k,C). Therefore, the number of positive solutions of (|10j) - (|lip is the same 
as the number of positive solutions of (|12p and (|14p . But v is uniquely determined by u in (|12p . which 
further simplifies the problem to one equation (I14p with one unknown u. Based on this observation, we 
have: 

Theorem 1 For each positive numbers S^ f, 7, there exist e$ > and k € R^"^ 6 such that the following 
property holds. Pick any E^gf, Ff f such that 

F tot = E toth < ^s tot h, (15) 

then the system S(k, C) with C = (Ef f, F^, Sf f) has at least n + 1 (n) positive steady states when n is 
even (odd). 
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Proof. For each k, 7, S^oti let us define two functions R+ x M + — > R as follows: 

F K 'T' 5 tot ( £) n ) = ^ K -( £5 tot' e5 tot/T- s tot) ( u ) (16) 

and 

F K ' 7 ' 5 tot ( e , u ) = F K ' {eS tot> £S tot^' S to\) ( u ) (17) 



= 2^(n) " 7 ^ (U)) 

By Lemma Q] and the argument before this theorem, it is enough to show that there exist Eq > and k € 
R^_ n_6 such that for all e € (0, £0), the equation F^'^tot (e, u) = has at least n + l (n) positive solutions 
when n is even (odd). (Then, given S tot , 7, E tot , and F tot satisfying ([15]), we let e = E tot /S tot < e , 
and apply the result.) 

A straightforward computation shows that when e = 0, 

i^< 5 tot( , u) = S tot (urf(u) - 7 ^(n)) - 7 ^(«) + u^(u) 

,n+l 1 \ \ I 1 , ^tot 



= A • • • A n _m" +i + A • • • A n _ 2 ^1 + (1 - 7&-1) - l^-ij u n 

+ • • • + A • • • Ai_ 2 (l + (1 _ _ 7 A i _ 1 J «* + •■■ (18) 

+ (l + ^(l-7/5o)-7Ao)^-7, 

where the Aj's and K^s are defined as in ©, and 0i = &cat Acat -- ^he polynomial F K ' 7 ' S tot (0, u) 
is of degree n + l, so there are at most n+l positive roots. Notice that u = is not a root because 
_p K,7,5 tot (0, u) = —7 < 0, which also implies that when n is odd, there can not be n + 1 positive roots. 
Now fix any S^ Q t and 7. We will construct a vector k such that F K,,y ' tot (0, u) has n + l distinct positive 
roots when n is even. 

Let us pick any n + l positive real numbers u± < ■ ■ ■ < u n+ \, such that their product is 7, and assume 
that 

(u - ui) ■ ■ ■ (u- u n+ i) = u n+1 + a n u n H h a\u + a , (19) 

where ao = —7 < keeping in mind that a^'s are given. Our goal is to find a vector k € R^j™ -6 such that 
(|18p and (|19p are the same. For each i = 0, . . . , n — 1, we pick Aj = 1. Comparing the coefficients of n 4+1 
in ([IHl) and {19]), we have: 

^tot, 



K Mi 

7ol ' 



'1 + a /3j) = a i+ i - a - 1. (20) 



Let us pick Km > such that f" M ' (ai + i — ao — 1) — 1 < 0, then take 

§rr( a i+i - °o - 1) - 1 

# = -£2s > 

ao 
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in order to satisfy (|20p . From the given 

Ao, • • • , X n -i,KM , ■ ■ ■ , Kj\i n -i,Po, ■ ■ • , Ai-ij 

we will find a vector 

« =(fcon ) • • • ,&on„_i, ^oif ' "" ' fc off„_i> fc cat ' • • • i*cat n -i' 
^on , • • • , ^on„_i , l o Q l ff n -i ' ^cat > • • • ' ^cat„_i ) G 11 



n6n— 6 



such that Pi = &eat Aeat > * = 0, ... ,n — 1, and (jSJ) holds. This vector k will guarantee that F K ' 7 ' 5 tot (0, it) 
has n + 1 positive distinct roots. When re is odd, a similar construction will give a vector k such that 
_p K ' 7,5 tot (0, u) has n positive roots and one negative root. 

One construction of k (given Aj, Km^ fa,i = 0, . . . ,n — 1) is as follows. For each i = 0, . . . , re — 1, we 
start by defining: 

XiKMi 

consistently with the definitions in Q. Then, we take 

^on, = 1, Wi, = 1) 

and 

1 — a 

fe offi = a i K Mi, k cati = (1 - ai)KMi, k&U = j3 , ~ K Mi, l Q Q. = L Mi - Z catj , 
where «j G (0, 1) is chosen such that 

I — a- 

l ofi l = L ^h 3~ KMi > °" 

This k satisfies fa = k ca ^_/l ca ^.,i = 0, . . . ,re — 1, and ©. 

In order to apply the Implicit Function Theorem, we now view the functions defined by formulas in 
(USD and ([I7D as defined also for e < 0, i.e. as functions RxR + — ► R. It is easy to see that F K,7 ' 5 tot (e, u) 
is C 1 on R x R + because the polynomial under the square root sign in F K '^' tot (e, u) is never zero. On 

the other hand, since F' l ' 7 ' s tot(0 , u) is a polynomial in u with distinct roots, — du ^ (0, Uj) 7^ 0. By the 
Implicit Function Theorem, for each i = l,...,re + l, there exist open intervals E{ containing 0, and open 
intervals Ui containing m, and a differentiable function 

a.i : E,i -> Ui 

such that Oj(0) = Uj, F^'^tot (e, o>i(e)) = for all e G Ej, and the images a^E^'s are non-overlapping. If 
we take 

n+l 

S "' -00) 



(0,e ) := fl ^fl (°' +c 



then for any e £ (0,eo), we have {aj(e)} as re + 1 distinct positive roots of J^' 7 ' tot (e, «). The case when 
re is odd can be proved similarly. I 

The above theorem shows that when -Etot/-^tot is sufficiently small, it is always possible for the futile 
cycle to have re + 1 (re) steady states when re is even (odd), by choosing appropriate kinetic constants k. 
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We should notice that for arbitrary k, the derivative of J- at each positive root may become zero, which 
breaks down the perturbation argument. Here is an example to show that more conditions are needed: 
with 

n = 2, A = l, Ai = 3, 7 = 6, ft = ft = 1/12, K = 1/8, K 1 = l/2, S tot = 5, 
we have that 

F^'^tot (0, u) = 3u 3 - 12u 2 + 15u - 6 = 3(tt - l) 2 (u - 2) 

has a double root at u = 1. In this case, even for e = 0.01, there is only one positive root of i^ re ' 7 ' tot (e, u), 
see Figure [2j 




Figure 2: The plot of the function F K,7,i tot (0.01, it) on [0,3]. There is a unique positive real solution 
around u = 2.14, the double root u = 1 of F^' tot(0, u) bifurcates to two complex roots with non-zero 
imaginary parts. 



However, the following lemma provides a sufficient condition for - — g- (0, u) 7^ 0, for any positive 



u such that F K ' 7 ' 5 tot (0, u) = 0. 

Lemma 2 For each positive numbers S^ot'T' an d vector k € R^j™ -6 , if 

1 - iPj - 1 



S tot 



< 



n 



(21) 



holds for all j = 1, • • • , n — 1, then 
See Appendix for the proof. 



dF K '~<' S tQt 
du 



(0,u) ^0. 



Theorem 2 For each positive numbers Sfgf,^, and vector k £ satisfying condition (|21|) . there exists 

£\ > smc/i £/taf for any F^ ^ Ef ^ satisfying Ff ^ = E^/j < the number of positive steady 

states of system S(k,C) is greater or equal to the number of (positive) roots of F re ' 7 ' foi(0,u). 



Proof Suppose that F* 1 '^' tot (0, u) has m roots: tii, . . . ,u m . Applying Lemma [21 we have 



Qp K n,s tot 
du 



(0,u k ) ^ 0, k = 1, . . . ,m. 
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By the perturbation arguments as in Theorem [H we have that there exists E\ > such that F K ' 7, tot(e, u) 
has at least m roots for all < e < e±. I 

The above result depends heavily on a perturbation argument, which only works when E+ ^/F+ + is 
sufficiently small. In the next section, we will give an upper bound of the number of steady states with no 
restrictions on ^tot/^tot) and independent of k and C. 

4.2 Upper bound on the number of steady states 

Theorem 3 For each k,C, the system S(k, C) has at most 2n — 1 positive steady states. 



Proof. An alternative approach to solving (I10p - (|llj ) is to first eliminate v from (jlOl) instead of from (lllj) . 
i.e. 

-^tot /-^tot ~~ u _ A(u) 

V = uifKu) - (E tot /F tot )^(u) := B{u)' 1 } 

when wpKu) — (^tot/^tot)^ ( u ) ^ 0- Then, we substitute (f22|) into (fTT|) . and multiply by (mp1(u) — 
(^tot/^totMKO 2 t° get: 

P K ' C ( U ) := <pff<p£ ( " «) 2 + fog - S tot¥ £ + *tot«tf + F tot ^) f |tot _ \ _ |tot ^ 

\*tot J V^tot / V not 

- Stat f utf - §2t ^ =0 . (23) 



Therefore, if wpi(u) — (-E'tot/^tot) ( / 3 2 ( n ) 7^ 0' the number of positive solutions of (fT0|) - (fTT]) is no greater 
than the number of positive roots of P K ' C (u). 

In the special case when u(pf(u) — (-^tot/-^tot)V2 ( u ) = 0; by (fT0|) . we must have it = ^tot/^tot' ana - 
thus tpl (-Etot/^tot) = ^2 (-^tot/^tot)- Substituting into (fTTj) . we get a unique v defined as in (fT2|) with 
it = -E'tot/^tot- But notice that in this case u = ^tot/^tot ^ s a ^ so a r0 °t °f P K ' C (u), so also in this case 
the number of positive solutions to (|10p - (|lip is no greater than the number of positive roots of P K,c (u). 

It is easy to see that P K,c (u) is divisible by u. Consider the polynomial Q K,c (u) := P K ' C {u)/u of 
degree 2n + 1. We will first show that Q K,c (u) has no more than 2n positive roots, then we will prove by 
contradiction that 2n distinct positive roots can not be achieved. 

It is easy to see that the coefficient of u 2n+1 is 

(A • • • A n _i) 2 Q 

and the constant term is 



E 



tot 



> 0. 



-^tot-^Mo 

So the polynomial Q K,c (u) has at least one negative root, and thus has no more than 2n positive roots. 

Suppose that S(k,C) has cardinality 2n, then Q K,c (u) must have 2n distinct positive roots, and each of 
them has multiplicity one. Let us denote the roots as u\, . . . ,ui n in ascending order. We claim that none 
of them equals E tot /F tot . If so, we would have </>? (^totMot) = ^(^totMot); and ^tot/^tot would 
be a double root of Q K ' C (u), contradiction. 
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Since Q K ' C (0) > 0, Q K ' C (u) is positive on intervals 

J = (0,Mi),Ii = (U2,U3), . . . ,I n -l = (u 2n -2,U2n-l),In = (u 2n ,Oo), 

and negative on intervals 

Jl = (ui,U 2 ), ...,J n = (u 2 n-l,U 2 n)- 

As remarked earlier, tp% (E tot /F tot ) ^ <p% (E tot /F tot ), the polynomial Q K ' C (u) evaluated at E tot /F tot 
is negative, and therefore, -^tot/^tot belongs to one of the J intervals, say J s = (u2 S ~i,U2 S ), for some 
s G {l,...,n} . 

On the other hand, the denominator of v in (|22p . denoted as B(u), is a polynomial of degree n and 
divisible by u. If -B(tt) has no positive root, then it does not change sign on the positive axis of u. But v 
changes sign when u passes £^ t/-^tot> thus V2 S -i and V2 S have opposite signs, and one of (u2 S -i, V2 S -\) 
and (u2s,V2s) is not a solution to (|10p - (|lip . which contradicts the fact that both are in S(k,C). 

Otherwise, there exists a positive root u of B(u) such that there is no other positive root of B(u) 
between u and E^/F^. Plugging u into Q K,c (u), we see that Q K ' (u) is always positive, therefore, u 
belongs to one of the I intervals, say It = {u2t, U2t+i) for some t G {0, . . . , n}. There are two cases: 

L ^tot/ F tot < u - We have 

«2,_1 < E ioi /F ioi < u 2t < u. 

Notice that v changes sign when u passes -^tot/^tot> so * n e corresponding V2 S -i and V2t have different 
signs, and either (u2 S -i, V2 S -i) ^ S(k,C) or (u2t,V2t) ^ <S(k,C), contradiction. 

2 - ^tot/^tot > «• We have 

u < u 2t +i < E tot /F tot < u 2s - 

Since v changes sign when u passes -^tot /-^tot ' so ^ n e corresponding V2t+i and t>2 S have different 
signs, and either (tt2t+i,«2t+i) ^ S(k,C) or (ti2 S ,^2s) ^ S(k,C), contradiction. 

Therefore, S(k, C) has at most 2n — 1 steady states. I 
4.3 Fine-tuned upper bounds 

In the previous section, we have seen that any (u,v) £ S(k,C),u 7^ ^tot/^tot must satisfy (f2"2|) - (f23"j) . but 
not all solutions of (i22l) - (i23l) are elements in S(k,C). Suppose that (it, v) is a solution of (|22l) - (|23"|) . it is 
in S(k,C) if and only if u, v > 0. In some special cases, for example, when the enzyme is in excess, or the 
substrate is in excess, we could count the number of solutions of (I22p ~ (l23p which are not in S(n,C) to get 
a better upper bound. 

The following is a standard result on continuity of roots; see for instance Lemma A. 4.1 in [30]: 

Lemma 3 Let g(z) = z n + a\z n ~ x + • • • + a n be a polynomial of degree n and complex coefficients having 
distinct roots 

Ai, . . . , X q , 

with multiplicities 

n\ H \-n q = n, 
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respectively. Given any small enough 8 > there exists a e > so that, if 

h(z) = z n + b\z n ~ l + ■ ■ ■ + b n , \cii — bi\ < e for i = 1, . . . , n, 
then /i has precisely roots in Bs(Xi) for each « = 1, . . . , q. 

Theorem 4 For each 7 > and k € IR^j™" 6 sitc/i that (pf(j) 7^ ^2(7); anc ^ each Sf f > 0, there exists 
£2 > such that for all positive numbers £^ i>^fof satisfying Ff f = E^/j < e^S^/j, the system 
E(k,C) has at most n + 1 positive steady states. 

Proof. Let us define a function R + x C — > C as follows: 

Q K ^ s tot( £ ,u) = Q K ' (e5 tot> £S tot^< 5 tot)( u ), 

and a set iJ^'^tot (e) consisting of the roots of Q K '^' tot (e, it) which are not positive or the corresponding 
v's determined by u's as in (f22j) are not positive, Since Q K 'T' S tot( £ , u) is a polynomial of degree 2n + 1, 
if we can show that there exists £2 > such that for any e € (0, £2), Q K ' 7 ' tot (e, u) has at least n roots 
counting multiplicities that are in S K ' 7 ' 5 tot( e ), then we are done. 

In order to apply Lemma [3l we regard the function Q K '' 1 ' S tot as defined on 1R x C. At e = 0: 

Q^< s tot(0,u) = [^(7 - uf + - Stot^XM - 7 ^)(7 " «) " Stot(«^ " 7^) 2 ]/« 

= [<^2(7 " u) 2 + <f${wpi - 7^(7 - u) - S tot (p%(u<p1 - 7^(7 - u) - S tot (u(p1 - ^) 2 }/u 
= [¥>o(7 " " ¥>a) + Stot«(«^ " 7^) (<^2 " ¥>?)]/« 

= tea - + StotM - 7^2) - 7^S) 

= (^-^)F^ S tot(0,n) 
Let us denote the distinct roots of Q^' tot (0, u) ju as 

Ul,... ,U q , 



with multiplicities 

and the roots of cpf — as 

with multiplicities 



ri\ H h n g = 2n + 1, 

p < q, 



m% + • • • + m p = n, nj > m,;, for i = 1, . . . ,p. 
For each i = 1, . . . ,p, if Ui is real and positive, then there are two cases ^ 7 as iff (7) 7^ ¥?2(7)) : 

1. Uj > 7. We have 

Ui<pt(ui) - Jip^Ui) > liipiiui) - ipZ{ui)) = 0. 

2. Uj < 7. We have 

UiP?(«i) - 7^2^) < j(<Pi(ui) - <pZ(ui)) = 0. 
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In both cases, unp1(ui) — 7</?2 (ui) and 7 — Uj have opposite signs, i.e. 

{ui<Pi(ui) - 7V2(«»))(7 - u < 0. 
Let us pick 6 > small enough such that the following conditions hold: 

1. For all i = 1, . . . ,p, if is not real, then Bs(ui) has no intersection with the real axis. 

2. For all % = 1, . . . ,p, if Uj is real and positive, the following inequality holds for any real u G B$(ui): 

(«^(«)-7^(«))(7-«)<0. (24) 

3. For all i = 1, . . . , p, if Uj is real and negative, then B$(v,i) has no intersection with the imaginary 



axis. 



4. B s (uj) fl Bs(u k ) = for all 3 + k = 1, . . . , q. 

By Lemma O there exists £3 > such that for all e G (0, £3), the polynomial Q K ' 7 ' tot(e,u)/u has exactly 
nj roots in each Bg(uj),j = 1, ... ,q, denoted by Uj(e), k = 1, . . . ,rij. 

We pick one such e, and we claim that none of the roots in B$(v,i),i = 1, . . . ,p with the v defined as 
in (|22p will be an element in S. If so, we are done, since there are Yli n i — X^i m i = n such roots, of 
Q K ^ s tot( £ ,u) which are in £ K ' 7 ' S tot ( e ). 

For each i = 1, . . . ,p, there are two cases: 

1. Uj is not real. Then condition 1 guarantees that u^{e) is not real for each k = 1, ... , rtj, and thus is 
in £ K ' 7 ' S tot( e ). 

2. itj is real and positive. Pick any root u^{e) G B$(ui), k = 1, . . . , rtj, the corresponding ^(e) equals 

7-^ fc (e) < Q 

(u?(e)¥>f(u?(e)) - 7<p%(u>{(e))) 

followed from (JS}. So («f(e),^(e)) g 5(k,C), and u*(e) G 5 K ' 7 ' 5 tot( e ). 

3. Uj is real and negative. By condition 1 and 3, (e) is not positive for all k = 1, . . . , rtj. 

I 

The next theorem considers the case when enzyme is in excess: 

Theorem 5 For each 7 > 0, k G R^j™ -6 such that (pi(j) 7^ ^2(7), and each Ef f > 0, there exists £3 > 
sitc/i t/iai for all positive numbers Ff f, Sf f satisfying Ff ^ = Efof/'y > Sf f/ '(£37), the system T,(k,C) has 
at most one positive steady state. 

Proof. For each 7 > 0, k G M+"'~ 6 such that ^1(7) 7^ ^2(7)' an d each £^tot > ^' we define a function 
R + x C — > C as follows: 

Q K ^ £ tot( e , u ) =Q K -( s tot' i? tot/^tot)(n). 
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Let us define the set C K, ^ ,E tot (g) as the set of roots of Q K '' y ' E tot (e ; u) which are not positive or the 
corresponding u's determined by u's as in (|22p are not positive. If we can show that there exists £3 > 
such that for any e € (0,^3) there is at most one positive root of Q K,1,E tot(e,u) that is not in C K,J,E tot(e), 
we are done. 

In order to apply Lemma [31 we now view the function Q K ' 7,S tot as defined on K x C. At e = 0: 

Q^^tot (0, u) = (7 - u) ((7 - u) rftf + + + ^<f£\ Ml - 7^2)) /« 

:= ( 7 -u) iT'^tot^). 
Let us denote the distinct roots of Q K ' 7,£; tot(0, as 

ttl(= 7),«2, ...,ttg, 

with multiplicities 

ni H \- n q = 2n + 1, 

and it2, . . . , u q are the roots of R K ' 1,E tot (u) other than 7. 

Since ^1(7) 7^ ^2(7)) R Kn ' E tot (u) is not divisible by it — 7, and thus ni = 1. 
For each i = 2, . . . , q, we have 

(E E \ 

Vo(«i) + + -^V2(«i)J («*Vi(w<) - 7¥>2 W) ■ 

If Ui > 0, then ^(it*)^ (ui) and ^o( n «) "I — ^r^UiPi H — ^P^fO"*) are both positive. Since Ui(pf(iii) — 
ytpziui) and 7 — Ui are non zero, unpf(ui) — "ff2( u i) an d 1 ~ u % m ust have opposite signs, that is 

(unpi(ui) - 7 ( / :, 2 (""■*)) (7 - u i) < 0- 

Let us pick 5 > small enough such that the following conditions hold for all i = 2, . . . , q: 

1. If Ui is not real, then B§(ui) has no intersection with the real axis. 

2. If Ui is real and positive, then for any real u £ B$(ui), the following inequality holds: 

W(«)-7V2(«))(7-«) <0. (25) 

3. If itj is real and negative, then Bg(ui) has no intersection with the imaginary axis. 

4. B (5 (« i ) n 5 5 (u fc ) = for alH ^ k = 2, . . . , q. 

By Lemma there exists £3 > such that for all e G (0,63), the polynomial Q K,n/,E tot (e, «) has exactly 
rtj roots in each Bg(uj),j = 1, . . . , q, denoted by Uj(e), k = 1, . . . , rij. 

We pick one such e, and if we can show that all of the roots in Bg(ui),i = 2, . . . , q are in C K ' 7,s tot (g), 
then we are done, since the only roots that may not be in C K,7 ' S tot(e) are the roots in Bg(u%), and there 
is one root in Bg(ui). 

For each i = 2, . . . ,p, there are three cases: 
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1. Ui is not real. Then condition 1 guarantees that u\ (e) is not real for all k = 1, . . . ,nj. 

2. ui is real and positive. Pick any root uf(e), k = 1, . . . , rij, the corresponding ff (e) equals 

7 ~ (e) < Q 

«i(^f(«?(£)) " 7^2 (^i( £ )) 

So, u*(e) is in C K ' 7 ' £ tot ( £ ). 

3. Uj is real and negative. By conditions 1 and 3, u^{e) is not positive for all k = 1, . . . , n%. 



5 Conclusions and discussions 



Here we have set up a mathematical model for multisite phosphorylation-dephosphorylation cycles of size 
n, and studied the number of positive steady states based on this model. We reformulated the question 
of number of positive steady states to question of the number of positive roots of certain polynomials, 
through which we also applied perturbation techniques. Our theoretical results depend on the assumption 
of mass action kinetics and distributive sequential mechanism, which are customary in the study of multisite 
phosphorylation and dephosphorylation. 

An upper bound of 2n — 1 steady states is obtained for arbitrary parameter combinations. Biologically, 
when the substrate concentration greatly exceeds that of the enzyme, there are at most n + 1 (n) steady 
states if n is even (odd). And this upper bound can be achieved under proper kinetic conditions, see 
Theorem Q] for the construction. On the other extreme, when the enzyme is in excess, there is a unique 
steady state. 

As a special case of n = 2, which can be applied to a single level of MAPK cascades. Our results 
guarantees that there are no more than three steady states, consistent with numerical simulations in [17]. 

We notice that there is an apparent gap between the upper bound 2n — l and the upper bound of n + 1 
(n) if n is even (odd) when the substrate is in excess. If we think the ratio E^/F^ as a parameter e, 
then when e<1, there are at most n + 1 (n) steady states when n is even (odd), which coincides with the 
largest possible lower bound. When there is a unique steady state. If the number of steady states 

changes "continuously" with respect to e, then we do not expect the number of steady states to exceed 
n + 1 (n) if n is even (odd). So a natural conjecture would be that the number of steady states never 
exceed n + 1 under any conditions. 
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7 Appendix 

proof of Lemma\^ Recall that (dropping the n's in ipf , i = 0, 1, 2) 

F K ^ 5 tot( ,u) = + S tot (wpZ - - 7Pg- 
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So 



du 

Since F K ' 7 ' S tot (0, u) = 0, 



(0,u) =ip% + Stot(«¥>i " 7^)' " (7 - u)(<p$)'- 



that is, 

„ _ - _ StotMi -7¥>2) 

! a — K 
^0 

Therefore, 

^ 5+ ^ w w-^)'-(, rf - w5 ) w )') 

= + -^((1 + ^o^ + A Aiu 2 + • • • + A • • • \ n -iu n ) x 

^0 

' (1 - 7 /3o) + 2-^(1 - j/3i)u + • • • + n A o -" A "- 2 (i _ 1 p n _ 1 )u n ~ 1 



Km Km-l ^Af„_! 
(A + 2A Ain H h nA • • • A 71— ill I X 

1 /- -n_ Aq , a \-2 , i Aq • • • A n _2 



(1 - 7 /3o)^ + -ttMI - 7&K + ' • • + " " (1 - 7/3n-iK ) 



ra-1 



+ % E A ° • • • A -^ Eg + 1 - o ^y^ a - 7/%)tf 

n 



^° i=o i=o 



'n-1 



+ Stot E A o • • • A -^ ( E^ + 1 - *°'" — (1 - 7&V) 

i=0 \j=0 M J 

= — E A ° " ' Ai -!^ A ° " " " A «-i^ n + E A ° " " " X i-^ j ( 1 + 5 tot(i + 1 
i=0 \ i=0 V 

where the product Ao • • • A_i is defined to be 1 for the convenience of notation. 
Because of (f2Tjh 

1-7/3, 



-Stot 



(i + i - *) 



K Mj 



< 1, 



, 9F K ' 7 ' S tot / n -n . n 

so we nave - — g- (U,it) > L). 
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